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Abstract 

It is well-known that in addition to the standard LMA solution to solar anomaly, there 
is another solution called LMA-Dark which requires Non-Standard Interactions (NSI) with 
effective couplings as large as the Fermi coupling. Although this solution satisfies all the 
bounds from various neutrino oscillation observations and even provides a better fit to low 
energy solar neutrino spectrum, it is not as popular as the LMA solution mainly because 
no model compatible with the existing bounds has been so far constructed to give rise 
to this solution. We introduce a model that provides a foundation for such large NSI 
with strength and flavor structure required for the LMA-Dark solution. This model is 
based on a new 17(1/ gauge interaction with a gauge boson of mass ~ 10 MeV under 
which quarks as well as the second and third generations of leptons are charged. We show 
that observable effects can appear in the spectrum of supernova and high energy cosmic 
neutrinos. Our model predicts a new contribution to the muon magnetic dipole moment 
and new rare meson decay modes. 


1 e-mail address:yasaman@theory. ipm. ac . ir 



1 Introduction 


Three neutrino mass and mixing scheme has been established as the standard solution to the 
solar and atmospheric neutrino anomalies. In general, the implicit assumption is that the 
alternative solutions to these anomalies {e.g., the magnetic transition moment of neutrinos) 
are now ruled out and can at most provide only a subdominant effect yet to be resolved. 
However, the so-called LMA-Dark solution, which is based on relatively large non-standard 
neutral current interactions of neutrinos with matter, is an exception that defies this general 
assumption. 

Surprisingly, the LMA-Dark solution not only passes all the bounds from various neutrino 
oscillation experiments but is also considered as one of the solutions for the suppression of the 
upturn on low energy spectrum of solar neutrinos which is in a mild tension with the standard 
neutrino oscillation scenario [1], In other words, the compatibility of LMA-Dark solution with 
the neutrino data is even slightly better than the standard three neutrino oscillation scheme 
without Non-Standard Interaction (NSI). The reason for this solution staying in the shadow of 
the LMA solution is mostly theoretical. In fact, there is no shortage of beyond standard models 
that give rise to NSI of neutrinos with matter [2,3]. The ratio of new effects within beyond 
Standard Model (SM) scenarios on forward scattering amplitude of neutrinos on matter to the 
SM contribution is expected to be given by {gx m x 2 ) / Gf , where gx and mx are respectively 
the coupling and mass of the new particle whose exchange leads to an effective four-Fermi 
interaction between neutrinos and matter field. Taking mx mw to avoid bounds from 
direct collider searches of new particles and requiring gx to be relatively small to remain in the 
perturbative region, we find that the NSI effects are suppressed relative to the standard model 
effects. For the LMA-Dark solution, the NSI effects should be as large as the standard ones. 
Thus, in the framework of beyond standard models with heavy particles {mx m-w), there is 
no theoretical basis for the LMA-Dark solution. The purpose of this letter is to build a model 
or a class of models that provide such a basis. 

The model presented here is based on a new 17(1) gauge interaction with a gauge boson 
of mass mz' ~ few 10 MeV and gauge coupling of g Z ' ~ 10 -5 which couples to the second 
and third generations of leptons (but not to the electron) as well as to the quarks. This gauge 
interaction leads to effective four-Fermi interactions with form and coupling strength required 
for the LMA-Dark solution. We expect new observable effects on the muon magnetic dipole 
moment, supernova evolution, trident neutrino production, meson decays and the spectrum of 
high energy cosmic neutrinos [4,5]. As shown in [6], such a model can be tested by studying 
/i + A —> + A + Z' where A is a nucleus. 

The letter is organized as follows: In sec. 2, the LMA-Dark solution is briefly reviewed. In 
sec. 3, the model is presented. In sec. 4, the potential signals of the model and observational 
bounds are reviewed. A summary of results is presented in sec. 5. 


2 The LMA-Dark solution 

Within the standard model, the neutral current interactions are flavor diagonal and universal. 
Going beyond the standard model, the neutral current interactions can have a more general 
form. Integrating out the heavy states, the correction to the four-Fermi neutral current inter- 
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action of the neutrino pair with matter can be described by the following effective operator 

C NS i = -2V2f ) ( 1 ) 

where / is the matter held (u, d or e), P is the chirality projection matrix and is a 
dimensionless matrix describing the deviation from the standard model. For propagation of 
neutrinos in matter, only forward scattering is relevant so integrating out the intermediate 
states in the /-channel amplitude and using the effective four-Fermi interaction is perfectly 
justified even if the intermediate states are relatively light. Since the atoms in the medium in 
which neutrinos propagate (the Sun, the Earth and etc.) are non-relativistic, mainly the vector 
part of the interaction contributes to the effective potential of neutrinos in matter which can 
be described by the following combination: 




+ e 


fR 

a.(3 * 


The non-standard interactions have been invoked in [1] to explain the suppression of upturn in 
low energy solar neutrino spectrum. This new solution which is called the LMA-Dark solution 
requires a value for #12 larger than 7 t/ 4 and 


- 1.40 < e“ e - < -0.68 and - 1.44 < e d ee - e d p < -0.87 at 3a C.L. (2) 

The off-diagonal elements of e Q( g and — e TT in the LMA-Dark solution can vanish. In [7], 
we have shown that the intermediate baseline reactor experiments such as JUNO can test this 
solution by determining which octant 9 12 belongs to 1 . NSI effects in long baseline experiments 
have been explored in [9]. Notice that the diagonal couplings of NSI at the LMA-Dark solution 
are as large as the standard ones. It is remarkable that such large values are still compatible 
with all the neutrino oscillation data [10]. However, as shown in [11], the results of the CHARM 
and NuTeV experiments rule out a significant (but not all) part of the range shown in Eq. (2). 
Notice however that to extract these bounds, it was implicitly assumed that the scale of new 
physics behind the effective couplings in Eq. (1) was much higher than the typical energy 
exchange in these scattering experiments rendering the effective Lagrangian formalism valid. 
If, as in the model presented in the next section, the masses of new particles are below typical 
energies at scattering experiments, these bounds do not apply. 


3 The model 

In this section, we show that it is possible to build a class of models that give rise to NSI of form 
required for the LMA-Dark solution (e MM — e ee ~ e TT — e ee ~ 1) avoiding all the present bounds. 
The models are based on a spontaneously broken 1/(1) gauge theory with gauge boson Z' p of 
mass of 0(10) MeV that couples to the second and third generations of leptons as well as the 
first generation of quarks. We denote the new gauge symmetry by U{1)'. To avoid the strong 
bounds on the coupling to the first generation of leptons, we assume that the 1/(1)' charges of 
the Erst generation of leptons vanish. Although coupling to the charged leptons or the second 

1 Notice that the reactor experiments are not directly sensitive to neutral current NSI but they can probe 
charged current NSI [8]. LMA-Dark solution does not however involve charged current NSI. 
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and third generations of quarks is not required for the LMA-Dark solution, we however need 
coupling to these fermions to cancel anomalies. In summary, the following gauge couplings are 
the basis of models that lead to NSI required for the LMA-Dark solution 


y l + y Qi QiY‘Q\ + r ul u R ru R + v d j R rd n ] z; e c 

«e(At,r} 


(3) 


where L a and Q i are the doublets of leptons of flavor a and quarks of first generation, re¬ 
spectively. Effects of couplings in Eq. (3) on the forward scattering (or any scattering with 
q 2 -C m|,) can be described by the following effective four-Fermi interactions: 


f-'NSI ~ 


g ,2 Y L 


m,-. 


^ {YQiQi') fjtQl 3“ Yui'U'R'yfi'U'R T ^di^R'y. 


(4) 


ae{n,r} 


Notice that we have taken the same couplings for /i and r. The mixing between v e and v [L (v T ) 
in the neutrino mass matrix is not compatible with the U'{1) symmetry under which z/ ;i and 
v T are charged but v e is neutral. The mixing of v e with and v T can be generated only after 
spontaneous breaking of the U'{ 1) gauge symmetry which requires new scalar(s). We will return 
to this point later. 

Equating C^si in Eqs. (1) and (4), we find 


and 


Requiring ~ 1, we find 


r „ s' 2 y l [y Qi +y„j 

TT ,Hl m%, 2V2 G f 

(5) 

* j f/ 2 Y L (Y Ql +Y dl ) 

TT NJ ' m|, 2V2 G f ‘ 

(6) 

g'~Tx 10~ 5 m . Z ' T . 

10 MeV 

(7) 


In order to avoid the bounds from big bang nucleosynthesis, the mass of mz' should be of order 
of 10 MeV or larger [4], Taking e“4 ~ 1, we find that at neutrino scattering experiments with 
four-momentum transfer of q, the contribution of new gauge interaction to scattering amplitude 
relative to the standard amplitude is given by m 2 z ,/(q' 2 — m|,). Notice that we have used the 
signature rf u = Diag(l, —1, —1, —1) so the t-charmel propagator is given by 1 /(q 2 — m 2 z ,). To 
avoid significant deviation from the SM prediction at scattering experiments such as NuTeV, 
NOMAD and etc., m Z ' should be much smaller than the typical energy of these experiments 
which is of order of GeV 

10 MeV ~ rriz' 1 GeV. 


Similar consideration holds valid for scattering of the muon beam (or hypothetical tau beam) 
on quarks. 

Up to now, we have only studied the couplings that give rise to NSI relevant for neutrino 
propagation. Let us now discuss the U( 1)' gauge couplings of all SM fermions and the possibility 
of the anomaly cancelation. Notice that Yl, (Yq 1 + Y U] ) and (Yq 1 + Y di ) should be all of the 
same sign to obtain positive value for effective couplings and defined in Eqs (5) and 
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( 6 ) as required for the LMA-Dark solution. As is well-known in the presence of right-handed 
neutrinos, the anomalies of £ e , C T and B /3 (where C Q is the lepton number of flavor a 
and B is the Baryon number) are equal. Thus, the combinations such as B — C, — C T or 
jC^+C t — (2/3)£> are anomaly free and can be gauged. In fact, there is rich literature on the B—C 
and — C T gauge theories. However, these models do not reproduce the flavor structure of 

and e d a o matrix as required by the LMA-Dark solution because while B — C predicts opposite 
sign for the [/(l)' charges of leptons and quarks, C M — C T predicts null coupling to quarks. 
Notice that the LMA-Dark solution differentiates between first generation of the leptons and 
the rest, ~ e'p? ^ e/f so we should assign different U(l)' charges to these leptons. Similarly 
to the lepton sector, let us define B, as quark number of generation i. Combinations of form 

Bfj, + B r + B\ — ciB -2 — (3 — a)B 3 ( 8 ) 

(with an arbitrary a) 2 * * * are also anomaly free. Gauging such a symmetry leads to 

Yq i = Y U1 = Y dl = 1/3 , Yq 2 = Y U2 = Y d2 = -a/3 and Yq 3 = Y U3 = Y d . 3 = -1 + a/3 (9) 

and 

n = Y Le = 0 and Y„ = Y r = Y w = Y Lt = 1 , (10) 

where Y e , Y tJ and Y r are U( 1)' charges of the corresponding right-handed leptons. Notice that 
with this assignment, the left-handed and right-handed fermions of each generation have the 
same U( 1)' charge so the diagonal elements of mass matrices do not break the U{1)' symmetry 
but the off-diagonal elements of the mass matrices (except the r/i element of neutrino mass 
matrix) break U( 1)'. As a result, to obtain mixing between different generations, new scalar(s) 
with proper [/(l)' charges should be introduced to couple to the fermions. The off-diagonal 
elements can be obtained after spontaneous breaking of electroweak and U( l ) 7 symmetries. 

Let us first discuss neutrino masses. The Dirac neutrino masses can be obtained from the 
following Yukawa couplings 

XiNiH T cL e + A 2 N 2 H T cL^ + X 3 N 3 H t cL t + \ a N 2 H t cL t + A 5 N 3 H T cL^ + H.c. (11) 

where Ni like L e is neutral under U(l)' but N 2 and N 3 are charged: Y N2 = Y N3 = Y L = Y Lt = 1. 
c is two by two antisymmetric matrix with c 12 = — c 2 \ = 1. Notice that by changing basis, either 
A 4 or A 5 can be set to zero but in general one of them remains nonzero and mixes the /i and r 
flavors with each other. However, to mix and v T with u e , we need extra scalars to break U (1)'. 
If we take only Dirac mass term, we shall encounter two problems: 1) Smallness of neutrino 
masses cannot be explained. 2) For Dirac neutrinos, we expect rn^. ~ m v . Considering that 
N 2 and N 3 are charged under U(l)' they can be produced via U( 1)' gauge interaction in the 
early universe and contribute to the relativistic degrees of freedom (A N e ff) on which there are 
strong bounds. To solve both problems, we can invoke the seesaw mechanism. Majorana mass 
terms for N 2 and N 3 require a new scalar (S/) with U( 1)' charge Ys 1 = —2. To mix v e with the 

2 More general combination of + cC T + bB\ — aB 2 — (1 + c + b — a)B 3 can also serve this purpose as long 

as arbitrary real numbers c and b satisfy c ~ 1 and b > 0. We however restrict ourselves to the less general 

combination shown in Eq. (8) to avoid unnecessary complexity. The case of c = —1 and b = —a has been 

proposed and explored in [12] to explain the LHC flavor anomalies. 
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rest we introduce another scalar S 2 with a U( 1)' charge of Y$ 2 = —1. We can then have the 
following terms in the potential 

A'hNfcNi + S 1 {A 2 N^cN 2 + A 3 N^cN 3 + A 23 NjcN 3 ) + S 2 (B 2 N?cN 2 + B 3 N?cN 3 ) + H.c. 


where Ai and Bj are dimensionless couplings. The mass of the Z' boson receives a contribution 
from the Vacuum Expectation Value (VEV) of Si and S 2 : .r/ 2 (V s ? (Si) 2 + Vj 2 (S 2 ) 2 ). However, 
since Si and S 2 are electroweak singlets, their VEVs cannot contribute to the masses of the W 
and Z gauge bosons. 

Let us now discuss mixing in the quark sector considering U{1)' charges shown in Eq. (9). 
To minimize the number of required scalars, we can choose the value of a such that the U{1)' 
charges of two generations of quarks become the same. It is straightforward to confirm that 
taking a = 3/2, a = 4 or a = —1 serves this purpose. Let us suppose that the U{1)' hypercharges 
of the i and j generations are the same but different from that of the k generation. Going to the 
mass basis, the ij element of the Z' coupling will be zero but the ik and jk elements will not 
vanish and can lead to flavor changing neutral currents. To avoid too large effects on the K — K 
and D — D mixing, we can set a = —1, rendering U(l)' charges of first and second generations 
the same. Thus, the Cabibbo mixing between Erst and second generations does not require 
breaking U( 1)'. To mix the third generation of quarks with the rest, we can introduce a second 
Higgs doublet H' with the same electroweak quantum numbers as those of the SM Higgs but 
nonzero U(l)' charge given by Y H > = Yq 12 — Y d3 = Yq 12 — Y u 3 = 5/3. With this new scalar, we 
can write Yukawa couplings of form H^cQ 3 u^ 2 and d 3 H^Q 12 - As a result, the VEV of H' will 
lead to the tc, tu as well as the bd and bs elements of the corresponding mass matrices but the 
ct, ut, db and sb elements of the same matrices vanish. The vanishing elements can be easily 
accommodated as we have the freedom to rotate the quarks with arbitrary unitary matrices. 
If the components of H' are not too heavy, they can be produced at the LHC and decay to 
quark pairs of different generations. Moreover, they can lead to new tree level flavor changing 
low energy effects. If H' is heavy enough (m^ 100 GeV), the bounds can be avoided but 

its VEV has to be smaller than (H) in the SM because, being electroweak doublet, it also 
contributes to the Z and W masses. This can be achieved by introducing another scalar, S, 
with U(iy charge of —5/3 and with the following potential 


- m 2 s \S\ 2 + \s\S\ 4 + msuH'SH^H’ + 


with 0 < m| -C rrijj, and m 2 EW < m 2 H ,. We will then have 



(H r ) = rn SH H' 


(H)(S) 

2 m 2 HI 


( 12 ) 


(13) 


Taking msHW , {H) -C rn h 1 and (S) ~ (H), we obtain (H') -C m^, itiew, (S)- Adding a 
quartic term for H' with a small coupling only slightly shifts (H'). Notice that there is a small 
mixing between S and H'° given by msHH' (H) /m 2 H , <Cl. A small mixing also appears between 
H' and H given by msHH' (S) /m 2 H , <C 1. The mass of Z' boson will be given by 

m Z ' = g' (Yl (S ,) 2 + Yl(S 2 ) 2 + y|(S ) 2 + Y 2 ,{H’) 2 ) in . (14) 
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Taking rriz' ~ 10 MeV and g' ~ 7 x 10 5 , we find that the VEV of electroweak singlet scalars 
should be of order of few hundred GeV. In particular, (Sf) rv_/ (S 2 ) r ''" / 100 GeV leads to right- 
handed neutrino masses of order of 100 GeV. Although, the masses are within the LHC reach, 
the g' coupling is too small to lead to a significant production of these particles at the LHC. 

If the H' components are not too heavy, they can be produced at the Run II of the LHC by 
electroweak interactions and lead to two jets of different flavors. 

4 Observational bounds and potential signals of the model 

In this section, we review the possible observational signals within this model and discuss the 
existing bounds. 

Meson decay: Taking Z' lighter than the pion, Hadronic decay modes for Z' will be forbidden 
and the tree level decays into and z/ T P r pairs will constitute the dominant decay mode. 

In this range of parameters, Z' can be produced in the meson decays and show up as missing 
energy. Some examples are K + —> fi + + + Z', K + —> e + + v e + Z', 7 r + —> g + + + Z' and 

7 r + —> e + + v e + Z'. The bound that has been obtained from the meson decay on the couplings 
of such a particle is of order of 10 ~ 3 [13] which is more than one order of magnitudes above the 
value we consider here (see Eq. (7)). Notice however that studying meson decays with higher 
accuracy can eventually provide a route to test these models. Since the branching ratios of 
these processes are proportional to g' 2 , to probe these effects precision on the branching ratios 
should be improved by a factor of ~ 200 . 

Mixing between Z' and 7 : Considering that the charged leptons as well as the quarks couple 
both to the photon and the Z' boson, a mixing will be obtained between ordinary photon and 
the Z' boson at one-loop level given by e ~ </e /(8 n 2 ). Such a mixing can affect the neutrino 
scattering on the electron. There are strong bounds from such interaction from Borexino 
experiment. The corresponding bound on the B — C models has been studied in [14]. Similarly 
to the C^ — C T model discussed in [4], the bound in our case has to be corrected by a factor 
of ( 1 / 0 . 66) 1 / 4 to account for the fact that only the z/^ and z/ r components of the solar neutrino 
flux couple to Z'. From the bound shown in Fig 8 of [14], we obtain ( g'ee ) 1//2 < 10~ 5 . Taking 
g' ~ 7 x 10 -5 , we find (g'ee) 1 / 2 ~ 2 x 10 6 which is well below this bound. The bound from 
GEMMA [14] does not apply for our model because z/ e does not couple to Z'. The beam dump 
experiments, which provide strong bounds on the B — C gauge theories are insensitive to this 
model because the produced Z' dominantly decays into vv pairs which, unlike the electron 
positron pair expected from the Z' decay within the B — C gauge models, leave no trace in 
detectors. In the parameter range of interest with rriz' ~ 10 MeV, there is no significant bound 
on the coupling of neutrinos to quarks [15]. 

Big bang nucleosynthesis: Decay of Z' into and v T v T pumps entropy into the neutrino 
gas which may appear as A N e ff. It is however shown in Ref. [4] that for > 10 MeV, the 
contribution to A N e ff is less than 0.1 and therefore negligible. 

Effects in Supernova: The effects of NSI in neutrino propagation in supernova and conse¬ 
quences for future supernova detection have been studied in literature [16]. Those results also 
apply to our model. Moreover Z l boson can also be produced inside the core and subsequently 
decay with ctz> ~ 10 _ 9 km(g'/7 x 10 _ 5 )~ 2 (T/10 MeV)(10 MeV/m ^/) 2 which is much smaller 
than the supernova radius. However, its production and subsequent decay can shorten the 
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mean free path of and v T and subsequently prolong the diffusion time [4], In case of future 
supernova detection, such effects can be resolved. A more detailed study of such effects is 
beyond the scope of the present paper. 

High energy cosmic neutrinos: The Z' boson can be resonantly produced via interactions of 
high energy cosmic neutrinos with background neutrinos which may cause a dip in the spectrum 
of cosmic neutrinos between 400 TeV to PeV [4,5]. In fact, the ICECUBE data shows hints for 
such a dip but further data is required to establish its presence. From the perspective of high 
energy cosmic neutrinos, this model is very similar to the — C r model studied in [4], As 
shown in [4], for the parameter range of our interest (i.e., mz> ~ 10 MeV and g' ~ 7 x 10~ 5 ) 
the optical depth can be larger than one and the dip can be therefore discerned. 

Muon magnetic dipole moment: The new 1/(1)' gauge symmetry leads to a correction to the 
muon magnetic dipole moment similar to the one-loop U(1)em correction. Up to corrections 
of 0(m|,/m 2 ) ~ 0 .01, we can write the contribution of the Z' loop as A(g — 2)^/2 = g' 2 /&n 2 
which for g' = 7 x 10~ 5 is below (but rather close to) the observational bound. This correction 
cannot account for the claimed discrepancy between SM prediction and measured value [17]. 

Signatures at the LHC: The components of the H l doublet can be produced at the LHC via 
electroweak interactions provided that they are not too heavy. Their decay will produce quark 
jets. They can also lead to flavor changing processes at tree level. However, we have introduced 
a mechanism to increase H' mass to arbitrarily large values to avoid bounds. 

Neutrino scattering experiments: We expect the effects of new physics at neutrino scattering 
experiments such as the NuTeV and CHARM experiments with energy momentum transfer q 
to be suppressed by a factor of m 2 z ,/\rn 2 z , — q 2 \ relative to the SM effects. Remember that if the 
scale of new physics responsible for the effective potential is much larger than q (as by default 
assumed in the literature), the correction to the cross section will be of order of main standard 
model. With this assumption, it was shown in [11] that a large part of LMA-Dark parameter 
space is ruled out by combining the results of NuTeV and CHARM scattering experiments. 
These bounds do not apply for our model because m 2 z ,/q 2 <C 1. 

Neutrino trident production: The new gauge interaction can contribute to the /i + /i _ produc¬ 
tion in the scattering of neutrino beams off nuclei which is known as neutrino trident production: 
v + A —> v + A + pA + pT where A is a nucleus. The CCFR collaboration [18] by studying the 
scattering of ~ 160 GeV neutrino beam off an iron target and the CHARM II collaboration [19] 
by studying the scattering of ~ 20 GeV neutrino beam off a glass target have extracted the 
cross section of the trident scattering process and found that the cross section is consistent 
with the SM prediction. From this observation, Ref. [20] has found a bound on g' which for 
rriz' ~ 10 MeV is around 9 x 10~ 4 . This bound is one order of magnitude above the values we 
are interested in. 

/i + A —>• p + A + Z', Z' —>■ vv: Z' can be produced by scattering of muon beam off nuclei 
and can then decay into a vv pair. Ref. [6] proposes using muon beam with energy of 150 GeV 
from CERN SPS to search for such a signal. It is shown that with 10 12 incident muons, values 
of g' as small as 10~ 5 can be probed which tests the present scenario. 
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5 Discussion and Summary 

We have presented a model leading to NSI required for the LMA-Dark solution. This model 
is based on a new U( 1)' gauge symmetry. To cancel anomalies, we have suggested gauging 

+ C T + B\ — aB 2 — (3 — a)B^ where C^ and C T are lepton numbers of flavor /i and r and 
Bi is the Baryon number of flavor i. The parameter a is an arbitrary real number. Taking 
rriz' ~ 10 MeV and g' — (6 — 8 ) x 10 _5 (m^//10 MeV), the values of and will be in 
the range required for the LMA-Dark solution (see Eq. (2)). Since the U{1)' charges of the g 
and r flavors are the same, the mixing between g and r flavors respects the gauge symmetry 
but mixing between v e and v tl (r ), or mixings between different generations of quarks break the 
U{ 1 )' symmetry. 

We invoke a seesaw mechanism in the neutrino sector and introduce two new scalars which 
are electroweak singlets but charged under the U(l)'. Their VEVs give Majorana masses to 
two right-handed neutrinos and mix v e with other flavors. In quark sector, we introduce a new 
scalar doublet H' with the same electroweak quantum numbers as those of the SM Higgs but 
charged under the U{ 1)' gauge symmetry. The new doublet has Yukawa couplings with quarks 
of different generations and its VEV mixes them. This new scalar doublet can contribute to 
flavor violation at tree level. To avoid the lower bounds from such processes and direct searches, 
H' has to be made heavy. We have suggested a mechanism to increase the mass of this new 
scalar while keeping its VEV small. The mass of Z' boson comes from VEVs of all these new 
scalars (see Eq. (14)). To obtain mz> ~ 10 MeV with g' ~ 7 x 10” 5 , these VEVs (or their 
maximum) should be of order of few 100 GeV. 

Obviously, this model, constructed to reproduce the LMA-Dark solution, has effects similar 
to what expected from the LMA-Dark solution on propagation of neutrinos in the Sun [1] 
and in supernova [16]. Moreover, as shown in [7], the solution can be tested by upcoming 
intermediate baseline reactor neutrino experiments through determining the octant of 612 . In 
addition to these effects, because of the presence of relatively light Z' boson with coupling to 
the standard model fermions with g' > 6 x 10~ 5 , we expect observable and testable effects in 
various experiments and setups: 

1. The correction to muon magnetic dipole moment is (g — 2)^/2 ~ 5 x 10” 11 which is too 
small to account for the claimed discrepancy. To probe such an effect, theoretical and 
observational uncertainties should be improved by an order of magnitude. 

2. The Z' boson can be resonantly produced via interactions of high energy cosmic neutrinos 
with the background relic neutrinos. This will create a dip in the spectrum of high energy 
neutrinos. The observed gap between 400 TeV to PeV might be due to such effects [4], 
After collecting more data by ICECUBE and by other neutrino telescopes, this prediction 
can be tested. 

3. The Z' gauge boson can be produced inside the supernova core. The produced Z' boson 
will immediately decay with a decay length much shorter than the supernova core radius. 

This can affect diffusion time scale of and T as well as the flavor composition and 
energy spectrum of supernova neutrinos. In the event of detection of supernova neutrinos, 
such predictions can be tested. 


4. The Z' boson can appear in meson decay as missing energy. The existing bounds constrain 
g' to be smaller than 10 -3 . To test the parameter range corresponding to the LMA-Dark 
solution, precision on the relevant processes such as [K + —> (/x + or e + ) + missing energy] 
should be improved by a factor of ~ 200. 

5. Z' can contribute to trident neutrino production: u + A —y u + A + /i + + ji~. From CCFR 
data, there is already a bound on g' [20]. To probe the range of parameter relevant for 
LMA-Dark solution, the bound should be improved by an order of magnitude. 

6 . The model can be tested by setup proposed in [6] which is based on studying fi + A —> 
fi + A + Z 1 , Z' —y vv with 150 GeV muon beam from SPS CERN. If the number of incident 
muons exceeds 10 12 entire parameter range relevant for the LMA-Dark solution can be 
tested. 

In summary, we have presented a UV complete model which gives rise to NSI required for 
the LMA-Dark solution based on a new U( 1) gauge interaction with a gauge boson of mass 
~ 10 MeV. Predictions of this model can be tested by future supernova neutrino observations, 
the spectrum of high energy cosmic neutrinos, trident neutrino production [20], study of rare 
meson decay, an order of magnitude improvement in measurement and calculation of (g — 2)^ 
and searching for p + A —> p + A + Z l via the setup described and proposed in [6] . 
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